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ABSTRACT 

Coherent  and  incoherent  ground-based  radar  measurements  of  the  winter  polar  cap  ionosphere  at 
Thule  and  Sondrestrom,  Greenland,  have  established  the  existence  of  "patches"  of  enhanced  ionization 
which  drift  across  the  polar  cap  in  an  antisunward,  noon-midnight  direction.  Associated  with  these 
patches  is  strong  radio  scintillation  activity  which  severely  disrupts  ground-to-satel 1 ite  communica¬ 
tion  systems  and  interferes  with  the  operation  of  space  surveillance  radar  at  high  latitudes. 

Several  recent  studies  have  shown  that  the  source  of  enhanced  ionization  is  the  sunlit  sub-cusp 
ionosphere  rather  than  production  by  precipitating  energetic  particles.  However,  the  question  of 
what  causes  the  "patchiness"  has  not  been  addressed.  We  study  this  problem  by  solving  the  time- 
dependent  plasma  continuity  equation  including  production  by  solar  ultraviolet  radiation,  loss 
through  charge  exchange,  and  transport  by  diffusion  and  convection  E  X  B  drifts.  Time  and  spatially 
varying,  horizontal  E  X  B  drift  patterns  are  imposed  and  subsequent  ionospheric  responses  are  calcu¬ 
lated  to  determine  how  enhanced  plasma  densities  in  the  dark  polar  cap  could  result  from  extended 
transit  of  relevant  flux  tubes  through  regions  of  significant  solar  production.  This  would  occur 
south  of  the  cusp  prior  to  convection  as  patches  across  the  polar  cap.  It  is  found  that  a  density 
enhancement  in  NMAX  from  7  X  104  to  5  X  10^  el/cm^  occurs  at  Thule  when  a  time-varying  convection 
pattern  is  included  in  the  simulation.  The  "patch  of  Ionization"  is  generated  when  an  initial 
convection  pattern  characterized  by  an  80  KV  crosstail  potential  and  a  12°  polar  cap  radius  is 
abruptly  changed  to  a  100  KV  crosstail  potential  and  a  15°  polar  cap  radius.  The  horizontal  extent 
of  the  patch  is  related  to  the  length  of  time  the  new  convection  pattern  remains  "turned-on". 

INTRODUCTION 

Substantial  evidence  exists,  both  observational  and  theoretical,  which  suggests  that  enhanced 
"patches"  of  ionization,  200  to  1000  km  in  horizontal  extent,  are  associated  with  a  southward 
directed  Interplantary  Magnetic  Field  (IMF)  orientation  (B2  <  0).  During  winter  hours 
(1000-2000  CGLT),  these  patches  are  observed  to  drift  in  an  anti-sunward  direction  over  Thule, 
Greenland  with  speeds  ranging  from  300  to  1000  m/sec.  Recent  observations  carried  out  simultan¬ 
eously  from  Thule  and  Sondrestrom,  Greenland  (Weber  et  al ,  1986)  have  traced  the  trajectory 
of  a  single  enhancement  for  over  3000  km.  Figure  1  presents  observed  fpF2  values  as  a  function  of 
Universal  Time  at  Thule  (86°  CGL)  for  the  period  4  Dec  to  10  Dec,  1983  (Buchau  et  al.,  1985).  The 
"patches"  of  ionization  are  evidenced  as  ionization  enhancements  normally  observed  between  1400  and 
2200  UT  and  may  be  broad  patches  as  on  5  Dec  or  "spikey"  patches  as  on  7  and  10  Dec. 

A  number  of  different  studies  have  demonstrated  that  the  plasma  within  the  enhanced  region  is 
produced  by  soiar  ultraviolet  radiation  in  flux  tubes  located  equatorward  of  the  subcusp  region 
before  convecting  over  Thule.  Buchau  et  al .  (1935)  have  presented  quantitative  arguments  which  demon¬ 
strate  that  the  observed  diurnal  variation  in  NMAX  (F2)  at  Thule  in  December  cannot  be  accounted  for 
by  local  production  due  to  precipitating  energetic  electrons  and  qualitatively  show  that  the  source 
of  plasma  is  equatorward  of  the  cusp  region. 

A  quantitative  ^ompa'-i 'on  of  calculated  and  observed  NMAX  (F2)  values  at  Thule  has  been  carried 
out  by  Anderson  et  al  (1986)  and  by  Buchau  et  al  (1985)  and  it  was  shown  that  indeed,  plasma  convec¬ 
ting  over  Thule  during  winter,  daytime  hours  (1000  to  2000  CGLT)  spends  about  3  to  5  hours  in  the 
sunlit  region  under  Bz  south  conditions.  Magnetic  flux  tubes  arriving  over  Thule  between  the  hours 
2000  and  1000  CGLT,  along  their  respective  trajectories  never  enter  the  sunlit  region  (see  fig  10, 
Buchau  et  al.,  1987,  IES  Proceedings)  and  ionization  production  is  due  entirely  to  precipitating 
energetic  electrons.  The  purpose  of  this  paper  is  to  present  a  possible  mechanism  for  creating  the 


649 


"patchiness"  in  the  plasma  observed  between  1400  and  2000  UT  and  to  answer  the  question  "wnere  are  the 
patches  formed?" 

THEORETICAL  CALCULATIONS 

Calculated  electron  density  profiles  as  a  function  of  latitude,  longitude  and  local  time  are 
found  by  numerically  solving  the  time  dependent  ion  (0+)  continuity  equation  given  by  expression  (1). 

3  Ni 

-  +  V  .  (NVi)  =  Pi-Li  (1) 

3  t 

Here  Ni  is  the  ion  density  (-Ne);  Vi,  the  ion  velocity;  Pi,  the  production  rate  and  Li,  the  ion  loss 
rate.  To  solve  equation  (1)  requires  transforming  the  independent  variables  from  a  spherical  r,  e, 
i>  coordinate  system  to  one  which  defines  directions  parallel  and  perpendicular  to  B.  After  some 
rearrangement,  equation  (1)  takes  the  form 

3  Ni 

-  +  Vi.  V  Nt  =  Pi-Li  -  V  .(NVi  1  )  -  Ni  .  V^  (2) 

3  t 

The  left  hand  side  of  expression  (2)  gives  the  time-rate-of-change  of  ion  density  in  a  reference 
frame  which  drifts  with  the  V^ convection  velocity  where  V  u=  E  X  B/g2.  A  detailed  description  of 
the  numerical  solution  of  expression  (2)  can  be  found  in  Anderson  (1973)  or  Moffett  (1979). 

We  assume  that  ionization  is  produced  by  solar  ultraviolet  radiation  of  wavelength  less  than 
91 1 A  and  neglect  the  contribution  by  precipitating  energetic  electrons.  Neglecting  particle 
precipitation  is  justified,  since  flux  tubes  convect  across  the  cusp  at  speeds  of  500  to  1000  m/sec-3 
(Heel  is,  19S4),  and  during  the  short  time  spent  within  the  high-flux  region  of  the  cusp,  the  expec¬ 
ted  contribution  to  the  resulting  electron  density  from  particle  precipitation  is  approximately  2.4 
X  104  cm"3  (Knudsen  et .  al.,  1977).  This  is  an  order  of  magnitude  less  than  the  densities  typically 
observed  within  the  patches  (3  X  105  cm-3,  or  foF2  =  5  MHz).  Loss  of  0+  ions  is  through  charge 
exchange  with  N?  and  0 2  where  the  loss  rate  coefficients  depend  on  the  relative  velocities  between 
ions  and  neutrals.  Plasma  transport  is  due  to  ambipolar  diffusion,  neutral  winds  and  convection 
E  X  B  drift  velocities.  Various  input  parameters  must  be  specified  in  order  to  solve  equation  2. 
These  models  are  now  described  briefly. 

1.  Neutral  atmospheric  densities  of  N^,  O2  and  0  and  the  neutral  temperature,  Tn,  as  a  function  of 
altitude,  latitude  and  local  time  are  obtained  from  the  MSIS  neutral  atmosphere  model. 

(Hedin  et  al.,  1977).  An  F10.7  cm  flux  value  of  120  is  chosen  to  represent  a  moderate  solar  active 
period  and  an  Ap  value  of  10  represents  quiet  conditions  in  January. 

2.  For  production  by  solar  ultraviolet  radiation,  the  photoionization  coefficient  at  the  top  of  the 
atmosphere,  P  „  is  chosen  to  be  5  X  10”'  sec  .  Plasma  diffusion  rates  are  similar  to  those  used  by 
Buchau  et  al  (1985)  and  the  loss  rates  due  to  charge  exchange  with  N2  and  O2  are  taken  from  Torr  and 
Torr  (1979).  The  appropriate  volume  rate  coefficients  are  dependent  on  effective  ion  temperatures 
which  depend  on  the  relative  motion  between  ions  and  neutrals. 

3.  Horizontal  E  1  I  convection  velocities  are  obtained  from  the  Heells  (1984)  model.  The  model 
specifies  the  geomagnetic  coordinates  (latitude  and  local  time)  where  the  plasma  drifts  to  after  a 
specified  At.  If  At  <0  then  the  model  specifies  where  the  plasma  came  from  At  seconds  earlier. 

Input  parameters  to  the  Heel  is  model  include  the  crosstail  potential,  the  width  of  the  "throat" 
region,  its  angular  displacement  relative  to  the  noon-midnight  direction  and  other  parameters  speci¬ 
fying  the  fall-off  of  potential  with  latitude  on  the  day  and  night  side. 

4.  For  this  particular  study,  neutral  wind  velocities  are  assumed  to  be  zero  since  they  are  of 
secondary  importance  In  actually  creating  the  "patches”.  In  reality,  however,  they  will  have  a 
moderating  effect  on  calculated  ion  and  electron  density  values. 

Plasma  which  crosses  over  Thule  at  a  specified  geomagnetic  local  time  is  assumed  to  originate 
six  hours  earlier  in  elapsed  time.  The  initial  location  in  geomagnetic  coordinates  is  found  by 
tracing  backwards  along  the  trajectory  established  (by  Heelis)  for  a  given  set  of  Input  parameters. 
After  each  tim<_-  otep  (300  sec),  the  appropriate  geographic  coordinates  and  solar  local  time  are 
calculated  corresponding  to  the  established  geomagnetic  coordinates  and  geomagnetic  local  time. 

Initial  locations  are  found  for  trajectories  which  pass  over  Thule  (six  hours  later)  at  hourly  Inter¬ 
vals  between  0800  and  2100  CGLT.  The  calculations  are  begun  at  these  initial  locations. 
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RESULTS  AND  DISCUSSION 


Two  convection  patterns  were  cnosen  for  this  study.  The  basic  differences  between  the  two  are 
character] zed  by  the  radius  of  the  polar  cap  region  and  the  magnitude  of  the  crosstail  potential. 

For  both,  the  throat  reyion  is  aligned  in  the  noon-midnight  direction.  One  pattern  possesses  a 
crosstail  potential  of  80KV  and  a  polar  cap  radius  of  12°  while  the  other  is  represented  by  a  100KV 
crosstail  potential  and  a  15°  polar  cap  radius. 

Figure  2  plots,  in  corrected  geomagnetic  latitude  and  lo  1  time,  the  trajectories  of  the  plasma 
which  cross  Tnule  at  1400  COLT  for  each  pattern.  When  the  8QKV  pattern  is  assumed  plasma  comes  from 
later  local  times  and  from  latitudes  of  70°  or  higher.  In  contrast,  when  the  100KV,  15°  radius 
Puttern  i  *  assumed,  plasma  originates  from  earlier  local  times  and  at  lower  latitudes  in  the  subcusp 
region.  The  time  between  points  on  each  trajectory  represent'  one  half  hour  elapsed  time. 

From  the  standpoint  of  plasma  production  by  solar  radiation,  however,  the  important  coordinate 
system  is  geographic  latitude  and  local  solar  time.  Figure  3  replots  the  trajectories  pictured  in 
Figure  2  in  this  coordinate  system.  The  solar  terminator  for  January  1  is  also  drawn  in  the  figure 
to  make  it  easier  to  see  which  parts  of  the  trajectory  enter  the  sunlit  region.  The  80K.V,  12°  radius 
trajectory  remains  essentially  in  darkness  the  entire  time,  barely  approaching  the  terminator  before 
crossing  over  Thule.  The  100KV,  15°  radius  path,  on  the  other  hand,  spends  more  than  five  hours  in 
the  sunlit  region  before  reaching  Thule.  This  is  due  to  the  fact  that  at  the  longitude  of  Thule  (and 
west)  a  lower  geographic  latitude  is  reached  for  a  given  geomagnetic  latitude.  For  longitudes  east 
of  Thule  at  later  local  times,  a  given  geomagnetic  latitude  corresponds  to  higher  and  higher  geo¬ 
graphic  latitudes. 

Ue  now  calulate  electron  density  profiles  as  a  function  time  by  solving  equation  (2)  numerically 
beginning  six  hours  before  the  plasma  reaches  Thule.  Figure  4  plots  calculated  NMAX  values  at  Thule 
as  a  function  of  magnetic  local  time  for  the  two  assumed  convection  patterns.  It  is  emphasized  that 
this  is  not  the  density  variation  along  the  path  traced  out  by  the  convecting  plasma,  but  is  rather 
the  calculated  density  when  the  plasma  reaches  Thule  at  these  various  magnetic  local  times  having 
traveled  along  different  trajectories. 

The  solid-lined  curve  gives  the  local  time  variation  of  NMAX  for  the  80KV,  12°  radius  convection 
pattern  and  demonstrates  that  calculated  densities  in  the  neighborhood  of  1  X  106  el/cm3  are  obtained 
between  1200  and  1900  magnetic  local  time.  Between  1900  and  2200  CGLT,  plasma  crossing  Thule  has 
entered  the  sunlit  region  for  a  brief  time.  Before  1100  and  after  2300  CGLT,  plasma  densities  all 
fall  below  2  X  104  el/cm3  indicating  that  the  convection  paths  remain  totally  in  darkness. 

If  ionization  production  by  precipitating  energetic  electrons  had  been  included,  NMAX  values  would 
have  been  maintained  above  5  X  104  el/cm3  (Buchau  et  al  ,  1985)  during  this  time. 

In  contrast,  calculated  NMAX  values  given  by  the  dashed-line  curve  for  the  100KV,  15°  radius 
pattern  reach  a  maximum  value  of  5  X  10^  el/cm3  between  1500  and  1900  CGLT  and  remain  above  1  X  10^ 
el/cm3  from  1100  to  2200  magnetic  local  time.  These  trajectories  remain  in  the  sunlit  region  a 
significant  amount  of  time  which  accounts  for  the  larger  density  values.  In  fact,  the  difference 
between  the  two  trajectories  is  most  sensitive  to  the  polar  cap  radius.  An  80KV,  13°  radius  trajec¬ 
tory  is  very  similar  to  the  12°  radius  path,  while  the  14°  radius  trajectory  is  very  similar  to  the 
15°  radius  pattern.  If  the  calculated  density  at  Thule  is  so  sensitive  to  the  polar  cap  radius,  what 
is  the  effect  on  NMAX  of  incorporating  a  time  varying  convection  pattern,  one  that  is  more  realistic 
than  the  steady-state  pattern?  Specifically,  we  want  to  answer  the  question  can  "patches"  of  en¬ 
hanced  ionization  be  formed  by  a  time-varying  convection  pattern? 
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We  have  calculated  values  of  NMAX  at  Thule  as  a  function  of  magnetic  local  time  from  1100  to 
2300  CGLT  where  the  convection  pattern  is  initially  assumed  to  be  the  80KV,  12°  radius  pattern.  At 
1700  UT  (1430  CGLT  at  Thule)  the  convection  model  is  suddenly  switched  to  the  100KV,  15°  radius 
pattern  and  maintained  for  two  hours  until  1900  UT  at  which  time  it  is  switched  back  to  the  original 
one.  In  the  second  case,  the  convection  pattern  is  changed  for  one  hour  from  1700  to  1800  UT. 
figure  5  displays  the  resulting  calculated  NMAX  values  at  Thule  for  the  two  cases.  When  the  100KV, 
15°  radius  pattern  is  maintained  for  two  hours,  the  change  in  NMAX  is  given  by  the  solid-lined  curve 
and  when  the  pattern  is  maintained  for  one  hour,  the  change  is  denoted  by  the  dotted-lined  curve. 

For  both  cases,  the  change  in  NMAX  beginning  at  1430  CGLT  is  dramatic.  Within  15  minutes  the  den¬ 
sity  increases  from  7  X  104  to  5  X  10$  el/cm3.  Keeping  the  100KV  connection  pattern  "on"  for  two 
hours  causes  the  ionization  patch  to  last  about  one  hour  at  Thule  before  decreasing  to  the  density 
value  characteri st ic  of  the  BOKv  pattern  at  1700  CGLT.  When  convection  is  changed  for  only  one 
hour,  a  "spike"  in  NMAX  occurs  between  1430  and  1530  magnetic  local  time. 

The  density  enhancement  at  Thule  responds  very  quickly  to  the  change  in  the  convection  pattern 
primarily  because  the  convection  velocity  is  so  high.  As  shown  in  Figure  3,  it  takes  plasma  only  1/2 
hour  to  travel  from  60°  geog.  latitude  to  Thule  (76°  geog.  lat.)  when  the  crosstail  potential  is 
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1 OOK V .  This  is  equivalent  to  a  convection  velocity  of  900  m/sec.  It  appears,  therefore,  that  a 
viable  candidate  for  creating  and  orginizing  polar  cap  patches  is  a  time-varying  convection  pattern. 

The  patch  is  formed  simply  by  allowing  plasma,  normally  found  in  the  sunlit,  sub-cusp  region,  access 
to  the  polar  cap  region  for  short  periods  of  time.  The  most  sensitive  convection  parameter  which 
permits  this  transition  to  occur  appears  to  be  the  polar  cap  radius.  The  angle  the  “throat"  makes 
with  the  noon-midnight  direction  does  not  seem  to  be  a  critical  parameter.  It  has  also  been  demon¬ 
strated  that  the  longer  the  convection  pattern  change  has  been  turned  "on",  the  larger  the  patch 
region  becomes  and  the  longer  it  takes  to  convect  over  a  given  location  in  the  polar  cap.  Ionization 
patches  which  are  more  “spikey"  in  nature  seem  to  represent  a  more  rapid  time-variation  in  the 
convection  pattern. 
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